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Abstract 

The Vliegenthart-Lekkerkerker relation for the second virial coefficient value at the critical tem- 
perature found in [G. A. Vliegenthart and H. N. W. Lekkerkerker, J. Chem. Phys. 112 5364 
(2000)] is discussed in connection with the scale invariant mean-field approach proposed in [V. L. 
Kulinskii and L. A. Bulavin, J. Chem. Phys. 133 134101 (2010)]. We study the case of the 
Mie-class potentials which is widely used in simulations of the phase equilibrium of the fluids. It is 
shown that due to the homogeneity property of the Mie-class potentials it is possible to connect 
the loci of the fluids with these model potentials in different dimensions. 
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I. INTRODUCTION 



Thermodynamic principle of the corresponding states (PCS) [fll |2j is one of the most 



vivid examples of the unifying nature of the scale invariance. From the microscopic point of 
view the PCS is based on simple scaling properties of the interaction potential which must 
be conformal ||. Actually, the microscopic interactions in real substances are more complex 
and do not conform with the conditions at which the PCS can be derived rigorously from the 
first principles of statistical mechanics. To expand the range of the applicability of the PCS 
its simple scaling form was extended to include more parameters which are connected with 
the most important properties of the interparticle interactions of the simple fluids (see e.g. 
||). The locus of the critical point and corresponding compressibility factor Z c = P c /n c T c , 
where P c , n c , T c are the critical pressure, density and the temperature correspondingly, is 
the core element of the PCS came originally from van der Waals. 

The interaction between nearest neighbors dominates in a dense, condensed phase of the 
molecular system. The liquid- gas critical point (CP) is adjacent to the region of condensed 
liquid phase. Therefore its locus is determined by the short scale properties of the interac- 
tions. They are commonly the distance which corresponds to the minimum of the potential 
and the depth the potential well e. But it is obvious that the very existence of the CP is 
due to the attraction at long distances || 0. For the potential $(r) of the general type 
its attractive and repulsive parts are independent. The short-ranged characteristics of the 
potential weakly depend on the long-range asymptotic behavior of the attractive part. For 
potentials of specific classes with simple analytical structure, e.g. the conformal potentials, 
it is possible to relate the short- and the long-ranged characteristics. In such case, as it was 
noted in 0, it is possible to construct the version of the PCS on the basis of a long-range 
interaction only. The coordinates of the critical point can be related with the long-range 
properties of the potential. For the potentials of the Lennard- Jones type such scaling ap- 
proach was proposed in ||. The generalized Lennard- Jones (LJ) potentials or Mie potentials 




q > m . 



(1) 
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serves as the example of conformal potentials. Here e is the depth of the potential well, the 
amplitude coefficient $o(<2S m ) is as following: 

MqM = - • (2) 



q — m \ m 

Below we use common dimensionless units for the temperature T — > T/e and the density 
n — > n o~ d , where d is the dimension. Due to advances in computer simulations there is a 
great body of results on the liquid-gas equilibrium for the systems with the potentials (P 



Obvious fundamental characteristic which incorporates both repulsive and attractive 
parts of the potential is the second virial coefficient B 2 (T) [12|: 

oo 

B 2 (T) = 2n j (l-e"*M/ T ) r 2 dr . (3) 
o 

In Ref. || it was demonstrated that for many three dimensional (3D) systems, where particles 
interact via spherically symmetrical pair interaction $(r), the relation: 

B 2 (T c )/v = -C , (4) 

holds, where C is some constant. It should be noted that C depends on the number of 
dimensions d too. As it follows from || the value of C weakly depends on the details of 
the potential as far as it belongs to some class of potential functions. In particular the au- 
thors of Ref. established that whereas the critical temperature T c drops considerably upon 
narrowing the range of attraction, the combination B2(T c )/vq remains practically constant 
0. According to for the generalized potentials $(r;6,g), C ~ 6.2 in three dimensions. 
These results put the question about the status of (^) as the mean of parametrization the 
classes of potentials in a sense more general than simple PCS. In such general statement 
the problem seems rather difficult. Nevertheless the possibility for its simplification may be 
found in reduction this problem for the continuous systems to their lattice analogues . 



In this paper we propose the derivation of the Vliegerthart-Lekkerkerker relation for the 
potentials of Mie-class (|l|) based on the global isomorphism between the Lennard- Jones 
fluids and the lattice gas model |§ . The potentials (0) have simple analytical structure and 
possess homogeneity dependence on the relevant parameters (m,q) (see Eq. (§)). 

The structure of the paper is organized as follows. In Section we discuss the dependence 
of the coordinates of the CP on the parameters of the potential based on the results of [EL [L4| . 



In Section |TTT| we consider the dependence of the Vliegerthart-Lekkerkerker parameter C on 
the characteristics of the Mie-class in d- dimensional case. The comparison with the known 
results is given and some predictions are made which can be tested in simulations. The 
discussion of the obtained results is in conclusive section. 



II. THE LOCUS OF THE CRITICAL POINT WITHIN THE GLOBAL ISOMOR- 
PHISM APPROACH 

It seems that there is no special symmetry property which governs the locus of the 
critical point in the continuum fluid. Because of the regular structure and the particle-hole 
symmetry of the configurations the lattice models possess additional symmetry properties 
which allow to calculate the locus of the critical point exactly without the direct calculation 
of the thermodynamic potential fl5| . For the continuum fluid the particle-hole symmetry is 
broken and the binodal is asymmetrical. The well-known rectilinear diameter law gives the 
representation of this asymmetry [pi pjfl. Recently the relation between the phase diagrams 



of the lattice gas and the simple liquid have been proposed in jy]] . It is based on the mapping 
between the fluid part of the phase diagram of the Lennard- Jones (LJ) fluids in the form of 
the projective transformation: 

n=n *TTIi' T=T *lfb (5) 

where x, t are the density and temperature variables of the lattice gas, n and T are the 
corresponding quantities for the continuum fluid. Here and T* are the parameters which 
will be defined below (see Eq. @). The temperature variable t is defined so that its value at 
the critical point (CP) is t c — 1. Here the parameter z is determined by the correspondence 
between the loci of critical points of fluid (n c , T c ) and the lattice gas x c — 1/2 , t c = 1: 

T, 



Z = ^-=r- (6) 



T*-T, 



c 



The inverse transformation has the form: 

n 1 T 

x= 1-t/t,' t = ~zf~~r- (7) 

The main assumption which governs the simple form of the transformation (|5]) is the validity 
of the rectilinear diameter law for the density in a broad temperature interval of liquid- 
gas equilibrium M. The inverse form (ffl) gives the procedure of the symmetrization of 



the binodal of the fluid in terms of the variable x which is combination of the density 
n and the temperature T of the continuum fluid. The idea that the difference between 
irregularity configuration for continuum fluids and regularity of configurations of lattice 
models is unimportant for order-disorder transitions was pronounced by K.S. Pitzer in [ITS] . 



It was stressed that the main difference is the shape of the holes in real or continuum liquid 
and that in the lattice gas. This is the reason why the behavior of continuum liquids beyond 
the fluctuational region is different from that of the lattice gas (Ising model). The global 
character of the transformation (||) shows that the particle-hole simplified picture is valid 
not only in near critical region but also far away from it. This rehabilitates the hole theory 
for expanded liquids [19] which gives the possibility to derive @ from the microscopic point 



of view. Yet the transformation Eq. (|D is not exact since its simple form is heavily based on 
the validity of the law of the rectilinear diameter. In the close vicinity of the critical point 



there are singular fluctuation corrections to the classical law of the rectilinear diameter [20 



Nevertheless, it is valid for the Lennard- Jones fluids in a broad temperature interval and is 



widely used to locate the critical point [[H|]. The numerical data of p2l |23[ are consistent 



with this law and therefore simple form of Eq. (Rl) is applicable for these systems. The 



results of |23J] for the mapping between Lennard- Jones fluids and the Ising model in 2 and 
3 dimensions support the validity of such approximation. 

From Eq. (f|) the loci of the CP for fluid and lattice gas are connected by simple relations: 



nc= 2(TT7)> Tc = T *— z > (8) 



where the parameter is given by: 



i?2 ( ) 

n * = T *BAf7)' (9) 



and T* is equal to the Boyle temperature in the vdW approximation T^ dW ^ = q/5 [^| |24 
Here b = 4t> is the fourfold of the eigenvolume of the particles v = 7rer 3 /6 and 



hoc 



a = -2tt / <£> attr (r)r 2 dr , (10) 



a 



where $ a ttr( r ) is the attractive part of the potential $(r). As has been shown in [14] it 



is possible to relate the parameter z of the transformation (Bf) with the exponent of the 



attractive part of the potential $ a ttr(?") — — r m . It is based on the relation 

dln(T c /r.) = 1 
d In ( ric/n* ) z ' 

which directly follows from Eq. (§) and the scaling relation 

d In ( T c /T* ) rf ln ( ^attr(^ c " 1/d )/T* 



(11) 



(12) 



d In n c / d In n c / 
between the critical temperature and that if $ a ttr( r ) — —r~ m . This relation generalizes 
commonly used reasoning that the critical temperature T c is of order of the magnitude of 
the potential well [[/J. Then in (^-dimensions z is determined as: 

z=±. (13) 
m 

According to Ref. M C ~ 6.3 in 3D for the potentials of the LJ type $(r;6, q),q > 7 as 
well as for $(r; m, 2m), m > 6. In Ref. |nj it was demonstrated that the results of Eq. (Bl) 



for the critical temperature are very close to those following from Eq. (f|) for the potentials 
with r -6 asymptotic of the attractive part. 

If the interaction does not have the hard core, there is an arbitrariness in the definition 
of the diameter of the particle a in Eq. fllP]) . Usually, it is chosen as the root of the equation 
$(cr) = 0. The locus of the CP determined by Eq. @ depends on the definition of soft core 
diameter a and changes with the change of this scale. Nevertheless, it should be noted that 
Eq. (H) does not depend on the spatial scale used to determine the size of the molecule like 
a. This size is introduced from the physical reasonings to distinguish between the short and 
the long range part of the potential (see e.g. |25|, |26|). 



Note that according to Eq. ([[]) the quantity $o is the homogeneous function of m and 
q. The change of the the scale for a leads to simple rescaling of the parameter T* because 
it is determined by the long-range power-like asymptotic behavior of the interaction. This 
points to the connection between (||) and Moreover, Eq. (||) in principle gives the value 
of the critical density too, while Vliegenthart-Lekkerkerker relation (|j) as it is formulated 
originally does not allow to estimate n c . Therefore it is interesting to clarify the interrelation 
between Eq. ([D and Eq. (H). We consider the results of Ref. || for the generalized Lennard- 
Jones potentials as the nontrivial extension of the PCS. The class of PCS consists of the 
family of potential functions given by Eq. ([!]). This inference is supported by the computer 
simulations [fLOl 01. 



Below we elucidate the nature of the Vliegenthart-Lekkerkerker relation (|4]) and its con- 
nection with the global isomorphism approach. We use the results of Ref. ]T1| to connect 
the value C with the parameter z. The latter represents the class of the thermodynamic 
similarity in accordance with the generalized principle of corresponding states JTJ1. 



III. THE NATURE OF THE VLIEGENTHART-LEKKERKERKER RELATION 
WITHIN THE GLOBAL ISOMORPHISM APPROACH 

The relation Eq. (^) looks quite unusual in a sense that it states the relation between the 
characteristic B 2 determined in low density region n — > with T c . The last is related with 
the moderate densities of liquid state. In the framework of the approach of Ref. [|17]] this is 
naturally explained since the Boyle temperature j ] ^ dW ' > — 7^ j s connected with the critical 
one T c according to Eq. (||). In previous works |14|, we show how the conception of the 



global isomorphism proposed in Ref. [|T7]] can be applied to formulate the scale invariant 
mean-field approach for calculation of the locus of the CP. 

We are not restricted by the 3D case and consider the case of general dimension d > 1 
so that the second virial coefficient is as following: 

00 

B 2 {T) = ^ J ( 1 - e -*M/* ) r^" 1 dr , (14) 


where Sd = 7r d / 2 /T ( ^ ) is the unit sphere in d- dimensional space. This gives us the possibility 
to establish the dependence of C on the relevant parameters in d- dimensional case. 

For the temperature we will use the common unit, the absolute value of the minimum e 
of the potential well. The spatial scale is commonly connected with the parameter a. But 
in accordance with the scale invariant nature of Eq. ([|) one can choose it in appropriate 
way. We use such freedom for the definition of the spatial scale which separate the intervals 
of distance where the repulsive and attractive forces dominate correspondingly. 

It is easy to derive the estimate for the parameter C using the high-temperature asymp- 
totic expansion for Eq. (|i~4]) (see e.g. p7| ): 



B 2 (T) = b(l- T*/T) + o ( T,/T ) . (15) 
Taking into account that in d dimensions b = 2 d ~ 1 Vq and using the relation (|J), from Eq. ([|) 
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we obtain: 

nd-l 

C^—. (16) 

For the LJ potential $(r; 6, 12) in 3D case for which z = 1/2 this gives C ~ 8, for 2D case 
2 = 1/3, C ~ 6 correspondingly. 

Below we improve the estimate (pj|) for C in a case of Mie potentials ([!]). For these 
potentials the left hand side of Eq. ([TJ]) transforms into: 

B 2 (T)/v = 2 d - 1 b , (17) 

where 

oo 

6o(T/$ 0) m/d, q/m) = J ( 1 - exp ( ^ (l - ) ) ) cfc , (18) 

o 

and x = r d and v = Sd/d (a/2) d is the d- dimensional volume. Notably, the parameter b 
depends on its arguments homogeneously. So if the critical temperature T c is determined by 
some constraint for B 2 (T) with the dimension d fixed, then it is the homogeneous function of 
the relevant parameters (m/d, q/m). We can say that the last represents the corresponding 
homogeneity class. The same homogeneous dependence is inherent to the parameter T* for 
the potentials of Mie-type ([]]): 

+ 00 q 

/( JL ) q-m 
$attr(r) r^ 1 dr = Kml . (19) 

(m/d — l)(q/m — 1) 

l 

Therefore, T* is the same for the potentials with the same ratio q/m. This allows to use sim- 
ple scaling considerations to connect the critical temperatures for the systems with different 
potentials belonging to the same homogeneity class. Besides, in view of the isomorphism 
with the lattice gas we can state that 

T c (d x )/T c (d 2 ) w d,/d 2 , (20) 

for the potentials of the same homogeneity class (m/d, q/m). Indeed, the expression @ 
assumes that the critical temperature of the isomorphic lattice model is set to unit t c = 1. So 
the correct comparison of the critical temperatures for the potentials in different dimensions 
demands the corresponding scaling since t c is proportional to the number of the nearest 
neighbors, which in its turn proportional to the dimension of the lattice d. Of course, 
the relation (pD[) is approximate since the fluctuations are neglected. Nevertheless it is 
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d 


2, $(r;6,12) 


3, $(r;9,18) 


4, $(r;12,24) 


T c 


0.5 


0.75 


1. 


rp(num) 
-t c 


0.515, @ 


0.73, § 


? 


n c 


0.353 


0.436 


0.55 


(num) 

rv c 


0.355, [28] 


0.354, § 


? 



TABLE I. The critical temperatures for the Mie potentials of the similarity class (m/d,q/m) = 
(3, 2). It includes $(r; 6, 12) in d = 2, $(r; 9, 18) in d = 3 and $(r; 12, 24) in d = 4. 



instructive to check the validity of Eq. (p0|) and the statement that the critical temperatures 
for the potentials belonging to the same similarity class can be connected via simple scaling 
relations. E.g. the potentials $(r; 6, 12) in d = 2 and $(r; 9, 18) in d = 3 fall into the same 
homogeneity class (2,2). Therefore 

T c {d = 3; 9, 18)/T c (d = 2; 6, 12) « 3/2 . 

Though the Vliegenthart-Lekkerkerker relation itself does not allow to calculate other pa- 
rameters of the critical point like the critical density and the pressure, they can be calculated 
using the methods developed in [Q . The value of the critical density is calculated in accor- 
dance with Eq. (P). Let us check (p0[) considering the potentials $(r; m, 2m) for d = 2,3, 4. 
For these potentials the parameter $o takes the constant value $o = 4 so that: 

B 2 (T) = 2 d - 1 b (T/A;m/d,2). (21) 



The results of calculations are presented in Tables |-|ITT| along with the known numerical 



estimates. We take the results for the critical temperature of the Lennard- Jones potential 
$(r; 6, 12) as known since they can be extracted from the corresponding lattice models (see 
also § H). 

Using relations (p|) and (0) for the case of the LJ potential $(r; 6, 12) we obtain: 
MT c /4; 2, 2) = SV lFl (-1; I; ?) r (?) - 3^F, g; |; |) T g) « 1.51 , (22) 



for d = 3. Here li^ is the Kummer confluent hypergeometric function |ftO|| . In accordance 
with Eq. (1171) this gives C ~ 6.04. In two-dimensional case for the LJ potential $(r; 6, 12) 



d 


2, $(r;4,8) 


3, $(r;6,12) 


4, $(r;8,16) 


T c 


0.89 


1.33 


1.78 


rp(num) 
J- c 


? 


1.313, 


? 


n c 


0.283 


0.322 


0.392 


(num) 


? 


0.316, [10] 


? 



TABLE II. The critical temperatures for the Mie potentials of the similarity class (m/d,q/m) = 
(2, 2). It includes $(r; 4, 8) in d = 2, $(r; 6, 12) in d = 3 and $(r; 8, 16) in d = 4. 



d 


2, $(r;3,6) 


3, $(r;9/2,9) 


4, $(r;6,12) 




1.6 


2.4 


3.2 


rplnum) 
J-c 


? 


? 


3.4, H 


n c 


0.333 


0.353 


0.41Q 


(num) 


? 


? 


0.34, || 



a In I we gave another value n c ~ 0.404. The difference is caused by different integration method and 
the precision for the calculation of multidimensional integral for B 3 in higher dimensions (d > 4). Here 
we obtained n c w 0.408 and rounded the result up to hundredth. 

TABLE III. The critical temperatures for the Mie potentials of the similarity class (m/d, q/m) = 
(3/2, 2). It includes 3, 6) in d = 2, $(r; 9/2, 9) in d = 3 and $(r; 6, 12) in d = 4. 

the value of the Vligerthart-Lekkerkerker constant C is determined by 

6 (T c /4; 3, 2) = f-i; 1; 2^ r (0 - 2 ^ Q; 5; 2 ) T « 3.81 , (23) 

which leads to C ~ 7.6. Finally, for the homogeneity class (3/2,2): 

MT c /4; 3/2, 2) . ( |) I ; 1) r (?) - V5 r g) ^ (1; 1 1)) » !,9 , 

(24) 

with the corresponding value C ~ 9.5. This corrects the high-temperature estimate fll6|). 
IV. CONCLUSIONS 



In this paper we derive the Vliegenthart-Lekkerkerker relation for the value of the second 
virial coefficient at the critical temperature [[| within the approach proposed in [14] for the 
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potentials of Mie-class. The very existence of such relation has appeared as the direct con- 
sequence of the global isomorphism between the fluids with the simple conformal potentials 
and the lattice models. We extended the applicability of this relation by considering the 
general case of d dimensions. The estimates for the loci of the critical points of the systems 
with $(r; m, 2m)-type potentials are obtained. These estimates can be checked in simula- 
tions and it would be interesting to investigate the relations between the critical parameters 
for the systems with the Mie potentials of the same homogeneity class but with different 
dimensions. This provides the test for the validity of the proposed approach. 

The application of the proposed approach to other types of the potentials with nonal- 
gebraic behavior, like hard spheres plus attractive Yukawa potential, Morse potential or 
square-well potential is connected with the derivation of the relation between the parameter 
z of the transformation (^) and the relevant integral characteristics of the potential. E.g., 
the relation ( |12|) can not be applied directly to the square well potential because of nonana- 
lytic dependence $(r). Obviously, these characteristics should be related with the symmetry 
properties of the potential under the scaling of its parameters, namely the effective range 
of interaction. The results of Ref. [[J show that the value of the Vliegenthart-Lekkerkerker 
parameter C differs essentially for Lennard- Jones like potentials $(r; 6, q) , q > 7 and square- 
well potential. In the last case there is strong dependence on the width of the potential well 
which governs the range of the interaction. The weak dependence of C on the exponent q of 
the repulsive part for $(r; 6, q) potentials in naturally explained within the global isomor- 
phism approach because of Eq. (p~3f ) and Eq. (§) . The square- well potential is the continuum 
version of the nearest neighbor interaction in the lattice model. So it is possible to connect 
the dependence of the critical parameters for the square-well potential and the correspond- 
ing results for the regular lattice gas model. These topics will be the subjects of further 
studies. 
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